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ABSTRACT 


This thesis discusses the modification of and additions 
to an existing Optimal Systems Control FORTRAN Program 
(OPTSYS) originally cktained from Professor Arthur E. Bryson 
of Stanford University. This program has been subsequently 
redesigned to run interactively on the IBM 3033 VM/CMS by 
lieutenant Commander John GG. Hoden, and additiors by 
Commander Harry A. Diel provide the user with a highly accu- 
rate graphic time response to a system designed using the 
Sevoxysk ELrogran. 

The addition of the FORTRAN program OPTGRAPH gives the 
user the capability tc obtain classical analysis (Pole-Zero 
Map, Foot Locus, Bode, Nyquist, and Nichols) of transfer 
functions, calculated for a system designed with the OPTISYSX 
program. The OPTGRAPH program uses high resolution precision 
plotting software to provide the user with highly accurate 
frequency response flcts. 

This series of programs permits users to rapidly carry 
out simulation, analysis, and design of Optimal Systems 
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I. INTRODUCTION 

The purpose of this theSis is to describe the modifica- 
tion and additions tc the existing FORTRAN program (OPTSYS) 
which is used in the study and application of Oftinal 
Systems Ccntrol theory. 

The OPtimal SYstems control program was originally 
developed by Hall [Ref. 1] to support his research in 
rotary-wing aircraft control systems. Later program modifi- 
cations were made by Walker [Ref. 2] and Liu [Ref. 3] of 
Stanford University, and are designated OPTSYS 4 and OPTSYS 
5 respectively. OPTSYS modifications made by Hoden [Ref. 4] 
were primarily devoted to creating a user-friendly interac- 
tive version (OPTS¥S3jNotf the OPTSYS 4 program: The latest 
modifications by Diel [Ref. 5] allowed the user to save 
matrices for subsequent runs of the OPTSYSX program, and 
formed data file sets for the time response program OPTCALC. 

The intent of this thesis work was to develop an inter- 
active program to plot the Pole-Zero map, Root-Locus, Bode, 
Nyguist, and Nichols responses to the open loop, closed loop 
noise, and compensator transfer functions of a State 
Variable Control System which has been developed using 
OPTSYS Program. Mincr modifications to the OPTSYS% Program 
were necessary to ccnstruct a data file sets for plotting 
the Pole-Zero Map, andin calculations for the Root-Locus, 
Bode, Nyguist, and Nichols plots. 

It is assumed that the user is familiar with the kasic 
concepts cf Control Theory and Optimal Systems Design. The 
synbol/naming conventions of Bryson [Ref. 6]. are used in 
the discussion of rogram operations and descriptions of 
problems using the OFISYS System. 


An overview of the OPTSYSX Program capabilities and of 
modifications to the existing programs is presented first. 
This is followed by a description of the program (OPTGRAPH), 
which was developed tc plot the Pole-Zero Map and to perform 
calculations for the Root-Locus, Bode, Nyquist, and Nichols 
Fiots. 

This work concludes with examples of various types of 
problems demonstrated in the interactive mode, including a 
copy of each terminal session with the final ‘results. 
Complete program listings for the OPTSYS EXEC progran, 
BeroYsk progran, and CPTGRAPH program are included in 


appendices A,B, and C respectively. 
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Ii. THE OPTSYSX GOMEUTER PROGRED 


A. BACKCROUND 


OFTSYSX is a double-precision, interactive FORTRAN 
program employing modern control theory analysis techniques. 
Its capabilities include the calculation of the cfen-loop 
eigensystem, and the stationary closed loop system; the 
synthesis of regulators and filters; along with pcwer spec- 
tral density, and modal distribution computations. The 
modifications introduced to the OPTSYSX program by this 
thesis work do net affect the program's original 


capabilities. 


Be OVERVIEW 


OFTSYSX is an extremely large and complex program with 
over 3000 lines of ccde. To use the program in its small 
versicn (dimensioned for a 32 xX 32 "F", "Gc", and ‘H" 
matrices), the user must extend his Virtual Machine (VM) 
memory capacity beyond the default VM memory capacity cf 720 
kilobytes to 1024 kilobytes. 

Any significant increase in the OPTSYSX program size, 
with the the resulting requirement for additional memory 
Capacity, would cause user difficulties. For this reason 
the task of obtaining the Pole-Zero maps, Root-Locus, and 
Bode, Nyquist, Nichols plots for the open loop, noise and 
compensator transfer functions was relegated to a separate 
program (OPTGRAPH). MThree data files (OPTGROL DATA, OPTGRNO 
DATA, and OPTGRCM DATA), containing the open loop, noise, 
and ccmpensator transfer functions, are formed by the OPTSYS 
program to be transfered to the OPTGRAPH progran. The 


OPTSYS EXEC program provides an interface between the 
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G@BTSYSX frogran, the time response programs (OPTCALC and 
OPTPLCT), and the OPTGRAPH progran. 


Camere rSYsk MODIFICATIONS 


CETSYSX program modifications consist of the addition of 
write statements, three flags, and a short routine to inter- 
pret user inputs. 

Write statements,to input data to the OPTGROL, CPTIGRNO, 
and CPTGRCM DATA files were added to the main frogran, 
Subroutine INNER, subroutine ZEROS, ands SuebOlimmne POLES - 
System information ccnsisting of number of states (Ns), 
number of controls (Nc), number of measurements (No), number 
of process noise scurces (Ng), type transfer function 
(ITFX), Markov Parameter {IE), and two flags to be used by 
the CFTSYS EXEC and the OPTGRAPH programs is inputed fron 
the main program. The write statements added to the subrou- 
tine INNER, and subroutine POLES input the poles’ for the 
open locp and noise transfer functions, and compensator 
transfer function respectively. The zeros, the numerator 
order, the gain, and the input and output numbers for all 
three transfer functions are obtained from the subroutine 
ZEROS. 

The first of the two flags, input with the system infcr- 
mation from the main frogram to the three data files, serves 
as a marker for the CPTSYS EXEC to locate the correct data 
dine for reading system inputs. The second flag, which is 
set by the routine te interpret user inputs, serves asa 
Signal for the OPTSYS EXEC to either load the OPTGRAPH 
program if any one of the three transfer functions was 
calculated by the OFTSYSX progran, and abort loading the 
OPTGRAPH program if none of the three transfer functions 
were calculated. Tiemeuleds flag {2TFX) is an existing 
G2Gowox Erogram flag. The flag (ITFX) is passed from the 


eZ 


subroutine TF to the subroutine ZEROS to idemtify the Cealen-— 
lation cf either an open loop, noise, Or compensator 
transfer function zeros. 

The short routine added to the main program of OPTSYSX 
reads the flags ITP, mit2, firs, wand hee The flag sent 
to the OPTGROL, OPTGRNO, or OPTGRCM DATA files is set to 
allow loading of transfer function data file to the CPTGRAPH 
program if the user selects either of the transfer function 
options number two (calculate poles, residues, and zeros - 
ITF1 = 1, &ITF2Z = AD oie ese) number three (calculate 
poles, and zeros - ITFI = 2, ITF2 =2, or, i283 — 2) xner 
calculating the open loop, noise, or compensator transfer 
functions. The selection of any other transfer function 
option will not calculate the transfer function zeros. The 
flag IRET is read to determine if the user desires to rerun 
a problem without exiting the OPTSYSX progran. Should the 
user elect to rerun a problem without exiting OPTSYSX (IRET 
= 1), the three data files are set back to the beginning 
(REWIND) to accept the new problem system data. 


De. OFTSYS EXEC MODIFICATIONS 


The OPTSYS EXEC was written by Diel [Ref. 53] to inter- 
face ketween the OFPTSYSX program andthe time resfonse 
programs OPTCALC and CPTPLOT. The modifications to the exec 
program were to add the option to select running the 
OPTGRAPH program to the exec menues and a short section to 
interface the OPTGRAPH program with the OPTSYSX progran. 

The added section to interface between the OPTSYSX and 
OPTGRAPH programs defined loader size (LDRTBLE), file defi- 
Nhitions (FILEDEF), and text libraries (TXTLIB) for the 
CPTGRAPH program. The first line of the three data files is 
read to ensure that there has been at least one of the 
three transfer functions calculated by the OPTSYSX prcegram 
prior to loading OPTGRAPH program. 


les 


The OFTSYS EXEC routine is also used to send _ flots to 
the VERSATEC plotter from the OPTPLOT program and the 
OPTGRAPH progran. 
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III. THE QPGRAPH PROGRAM 


A. PROGRAM OVERVIEW 


OPTGRAPH 1S an interactive FORTRAN program which uses 
the transfer function poles and zeros calculated by OPTSYSX 
to calculate, tabulate, and plot the Pole-Zero Map, Fode, 
Nyquist, Nichols, and Root-Locus responses for either the 
open loop, closed loop noise, or compensator transfer 


functions. 
1. Frogram Language 


CPFIGRAPH is programmed in FORTRAN following the 
conventions of the IBM System /360/370/ FORTRAN IV language. 
CPTGRAPH has been ccmpiled and run under both FORTRAN IV 
(G1) and FORTRAN H (extended) compilers on the IBM 3033. 


2- Graphics Package 


This program uses the Display Integrated Software 
System and Plotting System Software (DISSPLA) developed and 
distributed by the Integrated Software Systems Corporation 
(ISSCC) of San Diegoymealatornia. The DISSPLA package 1S a 
dibrary cf FORTRAN subroutines using FORTRAN IV conventions. 


3. Ilikrary Subroutines 


CPTGRAPH uses the International Mathematical and 
Statistical Library (IMSL) subroutine ZRPOLY™ tovecaleurae— 
the roots of a polynomial equation. ZRPOLY is a double 
precision FORTRAN subroutine to which is input a polynomial 
with real coefficients (double precision) in terms of 
decreasing powers of the polynomial. The subroutine finds 


the rects of the polynomial and returns a double precision 


ES 


complex array containing the roots. LREOLY is capable of 
accepting polynomials greater than 0 order and less than 101 


Cuder. 


CPTGRAPH haS one main program and 24 subroutines. 
The main program and its subroutines may be divided into 
four rasic categories: 
1) File Data Input 
2) Interactive Data Input 
3) Plot Setup and Sequencing 
4) Calculation 
A brief and general description of the program and 
its subroutines will follow in subseguent sections. 


B. GENERAL PROGRAM OPERATION 


CFTGRAPH was written to satisfy two specific objectives. 
The first is to create a program to be used aS a instruc- 
tmiondl tool for students taking controls related courses, 
and the second iS as method to assist in the analysis of 
actual large order ccntrol system problems. 

To accomplish these objectives the emphasis while 
writing this program was to make the program as user 
oriented as possible Ey attempting to eliminate ambiguities 
and providing features which would protect the user from 
inadvertent wrong entries. An effort was also made to mini- 
mize the Virtual Machine (VM) memory required to execute the 


program to facilitate the analysis of a large order systen. 


(eeoser Protecticn Features 


Three basic methods were used to protect the user 
from an abnormal program termination in the case of an inad- 


vertent wrong entry. The £1ESt waS to construct the 
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Subroutines to display a summery of the user inputs and 
allowing changes prior to leaving the subroutine. For the 
case where the user input is expected to fall within a given 
range, the input is examined and the user is either allowed 
to continue for an input within the range or issued an error 
or warning and reasked the question for an input outside the 
range. The final method is used by the three subroutines 
which read the user input from the screen. The subroutines 
RDINT, RDREAL, and RICHAR expect an integer input, real 
humber input, and a legiucal "Ves" OL Not wre speett yee In 
the event that a "null" line is inadvertently entered, the 
user is issued a warning and allowed another opportunity to 
enter the the correct input before abnormally terminating 
the progran. The entry of two "null" lines also offers the 
user the option to to exit the OPTGRAPH program at cther 
than normal program exit points. The subroutine RDCHAR also 
examines the input and issues a warning if the input is not 


either "yes" or "no", 
2. large Order Systen 


The large order control problem being considered for 
analysis by the OPTSYSX and the OPTGRAPH programs is the 
X-29A aircraft longitudinal axis backup mode systen. fail, 
system has a (98 X°9s) "PW matrix 7, "a (2 Xoo) “Hwee 
and a W998 xX “yee tater. 

To accommodate this large order system the OPTGRAPH 
program was dimensioned to accept a system with a maximum of 
99 states (Ns), 12 controls (Nc), 12 process noise sources 
(Ng), and 12 measurements (No). The Virtual Machine (VM) 
memory requirements for the program were minimized by the 
reuse of memory locations allocated for storage of calcu- 
lated data after that data has been either tabulated or 
plotted. This methcd of memory reuse has allowed the 
OPTGRAPH program to ke dimensioned to accommodate a large 
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order system and still operate with less than the 1024 
KILOBYTES VM memory reguired to run OPTSYSX in its small 


Vers 10n. 


Cw. SYSTEM/MNODEL DESCHIPTION 


The system equations used by the OPTSYSX program for the 
calculation of the transfer functions are of the state vari- 
able form. The system equations are: 
system model . 

xdot = [F]*x + [G]J*u + [GAN ]*w (ren) 
measurement equation 

2 = elt | D |e + “Vy 32) 
estimator equation 

RedQta- it Jaxer | Gatu + [ K |*(zZ = 1H |*xe) (Goce) 
open loor transfer function 

[iieis( fi > iF jl-2*[ 6] (3.4) 
closed lcop noise transfer function 

CHJ*is(l] - CF Ji-**[Gan] (3-5) 


compensator transfer function from measurement to input 


CCJ*CISEII-CFJ+(G J*0C J+CK J*C8 Ji-*) *[ 8] (3.6) 


where 


u = Control vector (Ns X 1) 
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W = White process noise vector (Ng X 1) 


x = State vector (Ns X 1) 

xdot = Derivative of the state vector (NS X 1) 
y = Output vector (No X 1) 

Z = measurement vector (No X 1) 


[Cc] = contro. Sarieiacei x, (Ne, Ns) 

[D}] = control feed-forward distribution matrix (Nc, Nc) 
of the control vector (u) 

[F] = open-loop dynamics matrix (Ns, Ns) 
(System matrix or Plant matrix) 

[G] = control distribution matrix (Ns,Nc) 

[GAM] = state disturbance distribution matrix (Ns,Nqg) 
[H1 = measurement distribution matrix (No, Ns) 


fIl1= identity matrix (Ns,Ns) 


De PROGRAM ORGANIZATION 
1. Main Program 


The main program presents a menu allowing the user 
to select for analysis either the oren loop transfer func- 
tion, the closed loop noise transfer function, the compen- 
Sator transfer functicn, or exit the OPTGRAPH program. The 
data file for the selected transfer function is examined by 
the Main program to ensure that the file contains the 
reguired transfer function data. Should the data file be 
incomplete the program gives the user a warning message 
indicating the nature of the missing data and allows the 
user te select another option. 


The main program menu: 
OFTIGRAPH 


DO YOU DESIRE TO ANALYZE: 
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1. OPEN LOOP TRANSFER FUNCTION 
2. NOISE TRANSFER FUNCTION 
3. COMPENSATOR TRANSFER FUNCTION 
G. EXIT OPTGRAPH 

ENTER OPTION NUMBER. 


Data is read from the selected file by the main 
program and formed into three arrays and seven integer vari- 
ables. The exchange of data between the main program and 
the subroutines, and between the subroutines is by the 
extensive use of labeled common statements. The data trans-~ 
fered through the subroutine call statements is limited 
primarily to flags, and constants. The program moves to its 
primary sequencing subroutine (GRAPH) upon the completicn of 


file data entry. 
2. fFrogram Sequencing 


The subroutine GRAPH first presents the user with a 
Summary cf the the transfer system received from the data 
fae . The user is then asked which transfer function he 
desires to to analyze using the conventions of the OPTSY5SX 
program describing the transfer functions by input number 
and output number. A menu is presented which allows the 


Mecmmeterseleet gGrapnical and/or tabular system response 


(Pole-Zero map, Root-locus, Bode, Nydiuiisteeeor Nichols) , 
select another transfer function, Or exit the OPTGRAPH 
ReOg i Alli. Upon completion of a graphical and/or tatkular 


response cption, or change of transfer function option, the 
program returns to the GRAPH menu to allow the user to 
choose another option or exit to the Main progran. 


Subroutine GRAPH data Summary and menu: 
OFIGRAPH 
THE OPEN LOOP TRANSFER SYSTEM OBTAINED FROM OPTSYS CONTAINS: 


1. NUMBER OF STATES = i 
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Z2- NUMBER OF CCNTROLS (INPUTS) = 1 
3. NUMBER OF MEASUREMENTS (OUTPUTS) = 1 
4G. MARKOV PARAMETER = 10.0**- 6 


CLEAR SCREEN TO CONTINUE 


OPTGRAFH 


AN OPEN LOOP POLE-ZERO, ROOT LOCUS, BODE, NYQUIST, AND/ 
CR NICHOIS PLOT IS 5u,EStREeeeon. 
INPUT # ? 


CUTPUT # ? 


OPTGRAFH 


OPEN LOOP TRANSFER FUNCTION 
INPUT # 1 
OUTEUT # 1 
DO YOU DESIRE GRAPHICAL RESPONSE AND/OR TABULAR DATA FOR: 
1. POLE-ZERO MAF 
2. ROOT-LOCUS 
3. BODE 
4. NYQUIST (POLAR PLOT) 
5. NICHOLS 
6. ANOTHER INPUT/OUTPUT COMBINATION (SAME TRANSFER 
FUNCTION) 
7. ANOTHER TYPE TRANSFER FUNCTION /EXIT OPTGRAPH 
ENTER OPTICN NUMBER. 


3- Graphic andfor Tabular System Response 


The selection of a graphical and/or tabular system 
response option for an Pole-Zero map, Root-Locus, Bode, open 
loop Nyquist, or ofen loop Nichols calls the subroutines 
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PZEno, si LOge BODE, NYCST, or NiCHOm respectively. While the 
response of these sukroutines is different, their crganiza- 
Lon 1S ‘similar. 

fie "SubprOuernes PZERO, RTLO, BODE, NYOST, and NICHOL 
primarily act as sequencing and plot setup subroutines. The 
interactive gquestioring of the user for graphic/tabular 
response information and the system response calculation has 
been relegated to other subroutines. The subroutines PZERO, 
.RTLO, BODE, NYQST, and NICHOL are divided into two major 
sections with the first being devoted to producing a graph- 
ical output and the second to producing a tabular output. 


The user is presented a subroutine menu offering the 


cption fcr graphic response, tabular response, or exiting 
the subroutine. The program moves to either the first 
section (graphic response), second section (je sushi iieie 


response) or back to the subroutine GRAPH menu depending on 
the option selection. The program returns to the subroutine 
menu upon completion cf the of either the graphic response, 
or takular response section, allowing the user to select 
another response opticn, or exit the subroutine. 

Features that are common to the five subroutines 
will ke described here and features that are unique to one 
of the subroutines will be described in the follcwing 
section under the aprprepriate response heading. 

For a graphic response, the user is given tke choice 
Meewe printer soptiens (TEK 618,.or. VERSATEC). With the 
selection of the TEK 618, the plot page size is defaulted to 
Mmimemrtcne= by S-5 tnches, and the plot will be presented on 
the TEK 618 screen. The VERSATEC option gives the user the 
option tc create up to a 21 inch by 21 inch plot (maximun 
POrminiemynkRoh! EC pPpLinter)) and Causes a DISSPLA METAFILE to 
be created. A scaling subroutine (PSCALE) scaies the plot 
heading, legend, etc. in proportion to the page width 
selected. The VERSATEC option offers the the advantage of a 


ZL 


high quality print” irom” the™ VERSATEOC Printer aot eee nc 
VERSATEC printer Opttenmin the OF15) sete theme re reno: 
a DISSPLA METAFILE also gives the user (after he exits 
optsys) the opportunity to use printers other than the TEK 
618 or VERSATEC. A major disadvantage with the VERSATEC 
cpticn is that the user must exit the OPTGRAPH program 
tefore he can print a graphic system response. For this 
reason, it is recommended that the system response first be 
pFlotted on the TEK 618 prior to selecting the VERSATEC 
option. 

For a tabular response,the user is presented a menu 
which gives him the choice of three devices (screen, 
printer, or disk) to send the tabular output. The user 1S 
asked if he desires to make any changes, after the tabular 
data has been calculated. If the answer is affirmative, he 
is presented a menu and allowed to make changes and rerun 
the tabular data routine. 

Examples of these features are included in the 
following sections with the interactive examples for the 


system responses. 


E. INTERACTIVE EXAMSFIES 


The open-loop, closed loop noise, and compensator 
transfer functions are are handled identically by the 
OPTGRAPH program, with the only difference being the program 
and graph headings identifying the transfer function tyfe. 
The frequency response programs for the Bode, Nyquist, and 
Nichols responses use identical setup and tabular data 
routines. 

A gocd overall representation of the OPTGRAPH program 
capabilities will be demonstrated by single interactive 
terminal session examples for the pole-zero, root locus, and 


bode responses using selected transfer functions. Tabular 


ae 


data examples will ke provided for the Pole-Zero, Root 
Locus, and Bode resfonses, and graphic response examples 
will bre provided the Fole-Zero, Root Locus, Bode, Nyquist, 
and Nichcls responses. 


F. CHARACTERISTIC EQUATION ANALYSIS 


The subroutines PZERO and RTLO are similar in their 
basic formats. The pole and zero locations for the selected 
transfer function infyut and output numbers are read into two 
Single dimension arrays. The extraneous zeros, calculated 
by the OPTSYSX progran, are eliminated by comparing the 
zeros to the Markov parameter sent with the system infcrma- 
Euan from the OPTSYSX progran. For a graphical response, 


the uSer is asked to define the plot limits in terms of 


X-coordinates (real axis) and Y-coordinates (imaginary 
mae S ) ls Points for the plot which fall outSide these limits 
will be ignored. System information (tranfer function 


input number, output number, and (DC) gain) is listed in the 


graph legend. 


1. Fole-Zero Mar 


The Subroutine PZERO converts the double precision 
numbers for pole and zero locations to single precision for 
compatibility with the with the Graphics Package and flots 
the poles and zeros within the plot limits defined by the 
user for the graphical response routine. For tabular data 
response the user is given system information (numerator 
order, denominator order, and transfer function (DC) gain) 
for the selected transfer function input number and outfut 


number, and the pole and zero locations. 
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The following example of a Simplified open loop transfer 
function for an airplane with an autopilot in the longitu- 
dinal mode was obtained from [Ref. 7,p.315]. 


The aircraft system is represented by: 


[(D] = [9] 


Oe eeee no ee Oo 
OReeeOs, "le Oem 
Je ce OS dca 
On) 16.-12. =3. |] 


[F] 


jy aor. 
[G] = | 0. 
i ae 
1-2. 


[GAM] = [0] 


[H) = — 1.90 oeno. ue 


a. Example cf Graphic Response (Terminal Session) 


BEGIN RECORDING OF TERMINAL SESSION 
Re T=020 170 s02 Vos 40: 17 
EXECUTION BEGINS... 


OPTGRAPH 


DO YOU DESERE BORA. 


ZS 


1. OPEN LOOP TRANSFER FUNCTION 
2. NOISE TRANSFER FUNCTION 
3. COMPENSATOR TRANSFER FUNCTION 
4. EXIT OPTGRAPH 

ENTER OPTION NUMBER. 

4 


CE TGHAPH 


Pit Oem he Gre mmanor ek SYSTEM OBTAINED FROM OPTSYs 


CONTAINS: 
1. NUMBER OF STATES = x 
Ze NUMBER OF CONTROLS (INPUTS) = q 
3. NOME EES Cr MEASUREMENTS \(OUTPUTS) = 1 


4. MARKOV PARAMETER = 10.0**- 6 


CLEAR SCREEN TO CONTINUE 


OPTGRAPH 


AN OPEN LOOP POLE-~ZERO, ROOT LOCUS, BODE, NYQUIST, 
Ripe OR NIGHOR eer cOtG>s DES@RED FOR: 
INP CL # = 


OUTPUT # ? 


OPTGRAPH 


OPEN LOOF TRANSFER FUNCTION 

INPCT # = 17 

OUTFOUT # 1 
DO YOU DESIRE GRAPHICAL RESPONSE AND/OR TABULAR 
DATA SOR. 

ie PeOLDE=ZERO MAP 
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2. ROOT-LOCUS 
3.  EODE 
4. NYQUIST (POLAR PLOT) 
5. NICHOLS 
6. ANOTHER INPUT/OUTPUT COMBINATION (SAME 
TRANSFER FUNCTION) 
7. ANOTHER TYPE TRANSFER FUNCTION /EXIT 
CPTGRAPH 
ENTER OPTICN NUMBER. 
: 


POLE-ZERC MAP 

OPEN LOOP TRANSFER FUNCTION 
DO YOU DESI EL: 

1. GRAPHICAL OUTPUT 

2. TABULAR DATA 

32) CU sUupROUTING 
ENTER OPTICN NUMBER. 
: 


POLEe=Z ERO eter 
OPEN LOOP TRANSFER FONCTION 
PLOTTER SELECTION AND PAGE SiZ2s 
(NOTE: PAGE IS 11.0 ENGHES SY S25 2NGHE Seer an 
SELECTIGN OFT EkG sp 
WHICH PLOTIER DO YCUuUCDESER.. 
Toe Eo nnS 
Zo VERSATEGC 
CHOCSE OPTION 1 OR 2 


PAGE SIZE SAMAXTMUMN = ZO eINGHE Ss 95 200 


INCHES) 
HEIGHT = 


JJ 


WIDTH = 
8 
PLOTTING LIMITS FOR GRAPH 
X AXIS (REAL AXIS) 
X MINIMUM = ? 


=) 
X MAXIMUM = ? 

5 

Y AXIS (IMAGINARY AXIS) 

YUEN = 2? 

=o 
Yon MM. 2 

5 


PLOTTING IIMITS FOR GRAPH 
X AXIS (REAL AXIS) 


X MINIMUM = alias 0.0 
X MAXIMUM = 5.00 
Y AXIS (IMAGINARY AXIS) 
Y MINIMUM = eee 
Y MAXIMUM = 5-00 


DO YOU DESIRE TO MAKE ANY CHANGES ? 
TYE Reyes eon UNO. 


n 
PLOT HEADING 


HOW MANY LINES OF HEADING DO YOU DESIRE ? 
( 3MAX ) 
ENTER NUMBER OF LINES. 
(ENTER "0" FOR NO HEADING) 


hetero Ofes2 CHARACTERS PER LINE IS 
ALLCWED 
Le 11S: 

Srnec himemw ten df auvcopi lot in 
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LING. 

the longitudinal mode 
LINE "Seite 

ref. 73; page 335 


PLOT HEALING 


LiN Eales ALRCRAFT WITH ANVAULOP Tie? ea 
LINE wees THE LONGITUDINAL MODE 
LINES ee REP so 7s PNGEaS 1 


DO YOU DESIRE TO MAKE ANY CHANGES ? 
TYPE "YES" Ofer. 
n 
>USING A PRE-ALLOCATED DATASET FOR UNIT FI1/F001. 
>USING A FRE-ALLOCATED DATASET FOR UNIT FT18F001. 


RK DISSPEA MEPAPILE HAS BEEN Create 
CLEAR SCREEN TO CONTINUE 


POLE-ZERO MAP 
OPEN LOOP TRANSFER FUNCTION 
DO YOU DESIRG®TO MAKE ANYMGHANGESe Te: 
1. PLOTTER /7 WAGES IZE 
2s GRAPH DPiuits 
3. HEADING 
4. NO ADDITIONAL CHANGES ~- PLOT POLE~ZERO 


MAP 
5. NO CHANGES = BXPTSrOrr-JErO PrP Ler PNG 
ROUTINE 
ENTER OPTION NUMBER. 


5 


POLE-ZERO MAP 
OPEN LOOP tRAN SEE reo Verbeet 
DO YOU DESERE- 
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ie 
2 
Bee 


GRAPHICAL OUTPUT 
TABULAR DATA 
CULT SUBROUTINE 


ENTER OPTICN NUMBER. 


3 


OPTGRAPH 


OPE LOOP TRANSFER FUNCTION 


oe 

Oye 
DO YOU DESIRE GRAPHICAL RESPONSE AND/OR TABULAR 
DATA FOR: 

1. FOLE-ZERO MAP 

2. ROOT-LOCUS 

3. EODE 

4. NYQUIST (POLAR PLOT) 

5. NICHOLS 

6. ANOTHER INPUT/OUTPUT COMBINATION (SAME 

TRANSFER FUNCTION) 
7. ANOTHER TYPE TRANSFER FUNCTION /EXIT 


CPTGRAPH 


ENTER OPTICN NUMBER, 


ii 


OPTGRAPH 


Pes YOU BEStRe ZOMANALYZEs 


US 
Bo 
35 
4 


Clee LOOP TRANSFER FUNCTION 
NOISE TRANSFER FUNCTION 
COMPENSATOR TRANSP Ener UNCTIION 
PA OP TGRAPH 


ENTER OPTICN NUMBER. 


4 
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R; T=2.80/9. Gye. 4515 
record off 
END RECORLCING OF TERMINAL SESSION 


The preceding example follows at the end of the 


chapter as figure Sai 
rE. Example of Tabular Data (Terminal Session) 


The entering and exiting option menus for the 
tabular data terminal session are identical to the graphic 
response terminal session option menus, and have been elimi- 
nated frcm the follcwing terminal session. 

BEGIN RECORDING OF TERMINALS SE sone! 
R; T=0.0 170.07 =o tee 
EXECUTION BEGINS... 


POLE-~ZERO MAP 

OPEN LCCE TRANSFER FUNCTION 
DOD OUD ES ERE 

1. GRAPHICAL OUTPUT 

Ze TABULAR DATA 

3. QUILT SUERCOUTE 
ENTER CPULCN SN Uae re 
2 


POLE-ZERO MAP 
OPEN LOCP TRANSFER FUNCTION 
DO YOU DESIRE TABULAR OUT oecee ler 
1. S@REBH 
Ze PRINTER 
3. DISK (OPGRAFH LISTING) 
ENTER CPIICN NUMBER. 
1 
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POLE-ZERO WMAP 
Srey LOOP TRANSFER FONCTION 


INEUT NUMBER = 1 
CUTEUT NUMBER = 1 
TRANSFER FUNCTICN (DC) GAIN = 0.1000D+01 
DENCMINATOR ORDER = 4 
POLE LOCATIONS 
REAL PART IMAGINARY PART 
0.0000D+00 0.0000D+00 
0.1000D+01 0.0000D+00 
-0.2000D+01 0.3464D+01 
-0.2000D+01 -0.3464D+01 
NUMERATOR ORDER = 1 
ZERO LOCATIONS 
REAL PART IMAGINARY PART 
-0.1000D+01 0.0000D+00 


Pome OU DESIRE TO CHANGE OUTPUT DEVICE? 
peee re ES AOR UNC. 


POLE-ZERO MAP 
OPEN LGOP TRANSFER FUNCTION 


wer Yeo DESIRE: 


1. GRAPHICAL OUTPUT 


2. TABULAR DATA 
Eee ours SUBROUTINE 


ENTER CPTICN NUMBER. 
3 


Pemem—teo0/5-25 19253200 


BEGCOEGd Off 
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END RECORDING OF TERMI s soso 


3. koot-Locus 


The subroutine RTLO creates a closed loop systen 
(with negative unity feed-back) from the selected transfer 
function input numbers and output numbers. Transier fune— 
tion roles and zeros are formed into two separate pclyno- 
mMials by the subroutine MAKPOL. The numerator polynomial 
(zeros) is multiplied by the gains, and the two polynomials 
are added to form a single polynomial. Complex roots of the 
of the pclynomial are calculated by the IMSL library sufrou- 
tine ZRECLY. The ccmplex roots are separated into real and 
imaginary components for either plotting or tabulating. 

The plotting routine first plots the transfer func- 
tion fole and zero lecations then plots the root locations 
as tkey are received from the subroutine ZRPOLY. For the 
fFlotting routine the gain interval specified by the user is 
divided into 2000 evenly spaced points for the subroutine 
ZRPOLY to calculate roots. The poles, zeros, and roots are 
plotted within the plct limits specified by the user. 

For the tabular data routine the user is given the 
cption te choose between 1 and 500 points for the sukrceutine 
ZRPOLY to calculate roots. The tabular data output provided 
for the user is system information (numerator order, denonmi- 
Nator order, and transfer function (DC) gain) for the 
selected transfer function input number and output number, 
the pele and zero locations, and the root locus gains and 
roots 

A positive unity feedback system may be analyzed by 


specifying a negative gain range vice a positive gain range. 


BS 


a.  Roct=Locus {Interactive Example) 


The following interactive example is of the ofen 
isepecradnsiter Luncticn for an aircraft with an autcpilot in 
the longitudinal mode described in Pole-Zero interactive 


secticn. 
a. Example of Graphic Response (Terminal Sessicn) 


muon RECORDING OF TERMINAL SESSION 
foe —O20170.02 20328337 
Ee SGCULION BEGINS... 


OPTGRAPH 


MemerCU DESTRE TO ANALYZE: 
1. OPEN LOOP TRANSFER FUNCTION 
Zee loros TRANSFER FUNCTION 
EeeCOUeENSATOR IRANSFER FUNCTION 
4. EXIT OPTGRAFE 

BitBR CPILCN NUMBER. 

1 


OPTGRAFH 


feeeCrhiy OOP TRANSFER SYSTEM OBTAINED FROM OPTSYS 


CONTAINS: 
1. NUMBER OF STATES = 4 
2. NUMBER OF CCNTROLS (INPUTS) = 1 
3. NUMBER OF MEASUREMENTS (OUTPUTS) = 1 
4. MARKOV PARAMETER = 10.0**- 6 


CrpAR SCR Reha hC CONTINUE 


OP VGh abn 
AN OPEN LOOP POLE-ZERO, ROOT LOCUS, BODE, NYQUIST, 
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AND/CHR NICHOLS PLOT “fee pes UR epee on. 
IN BUT 


OUTPUT # ? 


OPTGR APH 


OPEN LOCP TRANSFER FUNCTION 

INPUT # = 1 

CUTEUT # 1 
DO YOU DESIRE GRAPHICAL RESPONSE AND/OR TABULAR 
DATA FOR: 

1. POLE-ZERO MAP 

Zeos “ROCT=LOGUS 

Jon ODE 

4G. NYQUIST (POLAR PLOT) 

See NL GHO is 

6. ANOTHER INPUT/OUTPUT COMBINATION (SAME 

TRANSFER FUNCTION) 
va ANOTHER TYPE TRANSE ER FUNGRRONE na 


OPTGRAPH 
ENTER CPTION NUMBER. 
A 

ROOT -EOCUS 


OPEN LOCE TRANSFER FUNCTION 
TOn CU DEST aa. 
1. GRAPHICAL OUTPUT 
2. TABULAR DATA 
3. QULT SUSROUTINE 
ENTER CPIION NUMBER. 
| 


a 


ROOT-LOCUS 
OPEN LOOF TRANSFER FUNCTION 
FLOTTER SELECTION AND PAGE SIZE 
(NOTE: PAGE IS 11.0 INCHES BY 8.5 INCHES WITH 
SELECTICN OF TEK618) 


WHICH PLOTTER DO YOU DESIRE: 
1. TEK618 
2. VERSATEC 

CHCCSE OPTION 1 CR 2 


2 
PAGE SIZE (MAXIMUM = 21.0 INCHES BY 21.0 
INCHES) 
HEIGHT = 
6 
WIDTH = 
8 


PeOrranG LIMITS FOR CRAPH 
Menkes (REAL AXIS) 


X MINIMUM = ? 
aD 
X MAXIMUM = ? 
5 
Y AXIS (IMAGINARY AXIS) 
Veal Nr CO ee 
= 
Y MAXIMUM = ? 
5 


PLOPITING LIMITS FOR GRAPH 
X AXIS (REAL AX{MS) 
Reece Nt Wee = ed 
X MAXIMUM = a5 (00) 
Y AXIS (IMAGINARY AXIS) 


Io 


Y MINIMUM See iG 
Y MAXIMUM eu 
DO YCU DESIRE TO MAKE ANY “CHANGESS.: 
TYPE VYES” OR MSN e Le 


GAIN KFANCE FOR ROCi= Lecisee res 
MINIMUM GAIN = ? 


MAXIMUM GAIN = ? 
500 


GAIN RANGE FOR ROOGT= 1. CCUS es Ler 
MINIMUM GAIN 0.000D+00 
MAXIMUM GAIN 0.500D+03 
DO YOU DESIRE TO MAKE Aly cr ANGES 
TYPE MYES ORC ee 


EFLOT HEALING 
HOW MANY LINES OF HEADING DO YOU DESIRE ? 
( 3 MAX ) 
ENTER NUMBER OF LINES. 
(ENTER "0" FOR NO HEADING) 


A MAXIMUM OF 32 CHARACTERS PER LINE IS 
ALLOWED 
LINES WES: 
aircraft with an autoerior sin 
LUNE ees. 
the lengitudinal mode 
LINE Soc 
ref. 7; page 315 


FLOT HEALING 


3 


aN A 7s AIRCRAFT WITH AN AUTOPILOT IN 
Ne 2 .s Die LONGEiUDINAL MODE 
Sh a REE ee AGS 3 VS 
Dem OU DESTRE TO MAKE ANY CHANGES ? 
er eS en ANG, 
Nn 
PuotNnG A PRE-~ALLOCATED DATASET FOR UNIT FT1I7FOQ01. 
PUSING A PRE-ALLOCATEL DATASET FOR UNIT FTIT8FOO1. 


A DISSPLA METAFILE HAS BEEN CREATED 
CLEAR SCREEN TO CONTINUE 


ROOT-LOCUS 
OPEN LOOP TRANSFER FUNCTION 
WePYCU DESIRE TO MAKE ANY CHANGES TO: 
fee pee weOLlL ER 7sPAGE SIZE 
Ze GRAPH LIMITS 
3. HEADING 
4. GAIN RANGE 
See O ADDITIONAL CHANGES = PLOT ROGESsLOCUS 
Cees CHANGES. — EXIT ROOT-LOCUS PLOTTING 


Roe Oar N i 
EvlekR CPITCN NUMBER. 
6 
RCOT-LOCUS 


CPEN LOOP TRANSFER FUNCTION 
DO YCU DESIRE: 
1. GRAPHICAL OUTPUT 
2. TABULAR DATA 
o.oo SUBRCULINE 
ENTER CPIICN NUMBER. 
3 


ac 


OPTGR AEH 


OPEN LOOP TRANSFER FUNCTION 


INPUT # = 1 

Gule at o—n 
DO YCU DESIRE GRAPHICAL RESPONSE AND/OR TABULAR 
CATA FOR: 

1. POLE-ZERO MAE 

2. ROOT-LOCUS 

3. BODE 

4. NYQUIST (POLAR PLOT) 

5. NICHOLS 

6. ANOTHER INPUI/OUTPUT COMBINATION (SAME 

TRANSFER FUNCTION) 
7. ANOTHER TYPE TRANSFER FUNCTION /EXIT 


OPTGRAPH 


ENTER CPIICN NUMBER. 


7 


OPTGRAEH 


DO YOU DESIRE TO ANALYZE: 


OPEN LOOP TRANSFER FUNCTION 
NOISE TRANSFER FUNCTION 
COMPENSATOR TRANSFER FUNCTION 
EXD] OPITGRAPE 


ENTER CCPL CNR ER 


4 


Rs T=12- 367 85214) Joe eo 
TOECOEd OGL 
END KECOEDING OF TERMINAL SESSION 


a 


The preceding example follows at the end cf the 
enaoeecr aS figure 3.2... 


bk. Example of tabular data (Terminal Session) 


The entering and exiting option menus for the 
tabular data terginal session are identical to the graphic 
response terminal session option menus, and have been elini- 
nated from the fcllcwing terminal session. 

BEGIN KECORDING OF TERMINAL SESSION 
Pepe = Oe 70 .0220:42:59 
PRECUIICN BEGINS... 


ReEGI=LOCUsS 

Pen LOOF TRANSFER FUNCTION 
BOnLOU DESIRE: 

1. GRAPHICAL OUTPUT 

Ze 2 aewAR DATA 

nou Te SULROULINE 
Pie ek CPTICN NUMBER. 
2 


nROOT-LOCUS 
OPEN IOOP TRANSFER FUNCTION 
fee 2OU DESIRE TABULAR OUTPUT TO GO TO: 
Peer oN 
Zee eee NT ER 
3. DISK (OPGRAPH LISTING) 
Pee CPTLTCN NUNBER. 
1 


CeeneesanGh FOR ROOT=LCCUS PLOT 
MINIMUM GAIN = ? 


iH 
oJ 


MAXIMUM GAIN 
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GAIN RANGE FOR ROOT =~fCCiSeeLet 
MINIMUM GAIN = 0. 1000400 
MAXIMUM GAIN = 0.500D+01 
_2O YOU DESIRE TORI RE rAl ech iGi oe. 
TYPE “VYES “VOR Ge 


Nn 


HOW MANY PCINTS DO YCU WANT TO TABULATE? 


(500 IS THE MAXIMUM) 
ENTER, NUEBER Ol weCiy te. 


5 
ROOT-Leeus 
OPEN LOCF TRANSFER FONCTION 
INPUT NUMBER = 1 
CUTFUT NUMBER = 1 
TRANSFER FUNCTICN (DC) "GAIN = SOR CCCrs | 


OPEN ICCF TF POLES ANI ZEROS 


DENCMINATOR ORDER = U 
FOLE LOCATIONS 
REAL PART IMAGINARY PART 
0.0000D+00 0.00005+00 
0. 1O00D + On 0.0000D+00 
=J.20000 700 0.3464D+01 
-0.2000D+01 -0.3464D+01 
NUMERATOR ORDER = 1 
ZERO LOCATIONS 
REAL PART IMAGINARY PART 
-0.1000D+01 0.0000D+00 
GAIN = 0. 1000L5+01 
REAL PART IMAGINARY PART 
0.7074D-01 0.0000D+00 
0.8962D+00 0.0000D+00 
Ui oosDEon 0.3441D+01 


4 


-0.1983D+01 


-0.3441D+01 


GAIN = 0.18005+01 

REAL PART IMAGINARY PART 

0.1453D+00 0.0000D+00 

0.7948D+00 0.0000D+00 
-0.1970D+01 0.3422D+01 
-0.1970D+01 -0.3422D+01 
GAIN = 0.2600D+01 

REAL PART IMAGINARY PART 

0.2576D+00 0.0000D+00 

0.6552D+00 0.0000D+00 
-0.1956D+01 0.3402D+01 
-0.1956D+01 -0.3402D+01 
GAIN = 0.3400D+01 

REAL PART IMAGINARY PART 

0.4426D+00 0.1661D+00 

0.4426D+00 -0.1661D+00 
-0.1943D+01 0.3383D+01 
-0.1943D+01 -0.3383D+01 
GAIN = 0.4200D+01 

REAL PART IMAGINARY PART 

0.4285D+00 0.3096D+00 

0.4285D+00 -0.3096D+00 
-0.1929D+01 0.3363D+01 
-0.1929D+01 -0.3363D+01 


PemrecueURSIRE TO MAKE ANY CHANGES ? 


Peewee YES? OR "NOS. 


FOCT-LOECUS 


OPEN ICCF TRANSFER FUNCTION 
DOeYOU DESTRE TO MAKE ANY CHANGES TO: 


ie sOUlTPUT DEVICE 
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2. GAIN RANGE 
3. NUMBER OF POMS GAL Unies 
4. NO ADDITIONAL CHANGES ~- TABULATE DATA 
5S. EXIT ROOT=EPCCUS = IABULAr 2aia ROUSE 
ENTER CFIICN NUMBER. 
2) 


ROOT -LECus 

OPEN LOOF TRANSFER FONCTION 
[LO JYOUs DES Gre: 

1. GRAPHICAL OUTPUT 

2. TABULAR DATA 

3. OUT SUBROUTINE 
ENTER CPTICN NUMBER. 
S 


R; T=1.63/3.40 20:44:54 
FeCOEdwoOL: 
END RECOFDING OF TERMINAL SESSION 


G. FREQUENCY RESPONSE ANALYSIS 


The calculations for the frequency response of the Bode, 
Nyquist, and Nichols responses are identical. The only 
difference is in the graphical presentation of the calcu- 
lated data. The subroutines BODE, NYQST, and NICHOLS call 
the subroutine FREQ to calculate the transfer function 
magnitude and phase fcr the frequency range specified by the 
user. 

The subroutine FREQ eliminates the extraneous zeros by 
comparing the zeros to the Markov parameter as described in 
the Characteristic Equation Analysis section. The frequency 
response magnitude and phase are calculated using double 


precision complex libeass, functions following the 
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conventions of the IBM System /360/370/ FORTRAN IV language. 
The frequency response is calculated by determining the 
product of succesSive numerator zeros which have been 
combined with the frequency and divided by the preduct of 
successive denominatcr poles which also have been combined 
with the frequency. The magnitude for the complex result of 
the zeros, poles, and frequency calculations is determined 
using the FORTRAN dcuble precision library functions for 
finding the magnitude of a complex argument. The phase of 
the of the complex result is determined by separating the 
real and imaginary parts of the complex results, and finding 
the phase angle using the double precision FORTRAN library 
function arctangent. The angle is converted from radians to 
degrees and is corrected for proper gquadrant by checking the 
signs of the real and imaginary parts. The normal range for 
calculating phase angles by the subroutine is from +540 
degrees to -540 degrees. This range is considered adequate 
fOr most application, but should these limits be exceeded 
the flot will jump with a vertical line to either +360 
degrees or -360 degrees respectively before continuing with 
the plot. 

For a graphic response the user specified frequency 
range is divided into 500 equally spaced points (on a 
logrithmic scale) for which to calculate the magnitude and 
phase. The user is given the option to select between 1 and 


Boom points tome tabular data. 


1. Eode Response 


The graphic response section of the subroutine BODE 
creates two plots. The first plot is the response magnitude 
in decibels versus the frequency in radians per second. The 
Magnitude from the subroutine FREQ is converted to decibels 
and plotted with its corresponding frequency. The second 


plot is of the phase in degrees versus frequency in radians 
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fer second. The phase from the subroutine FREQ is plotted 
with its corresponding frequency. The frequency for both 
flots is plotted on a logrithmic scale. The scale range for 
both is set automatically to include the entire range of 
data to ke plotted by the subroutine MINMAX which scans the 
data for the maximum and minimum values. System information 
(tranfer function input number, output number, and (ie 
gain) is listed in the graph legend. 

The magnitude in decibels and the phase in degrees 
1s searched for "0" decibel and 180 degree crossover feints 
to calculate the phase and gain margins. 

The tabular response section gives the user the 
opticn to select the the number of frequency points’ to be 
calculated between 1 and 500 for the frequency range he has 
specified. For takular data response the user iS given 
system information (numerator order, denominator order, and 
transfer function (DC) gain) for the selected transfer func- 
tion input number and output number, and the frequency 
(cadians per second), magnitude (decibels), phase (degrees), 


and the real and imaginary frequency response parts. 


2. Eode Response (Interactive Example) 


The fcllcwing example of compensator transfer function with 
a filter and regulator was obtained from [Ref. 8,pp. 382 - 
384 }. 


The compensated system 1S represented by: 


[C] = { -254.1 -19.57 } 
L Oe 
Ceo 6, 1.0 j 

| O. -4.6 | 


[GAM] 


CH] 


[K] 


[QJ 


Example cf Graphic Response (Terminal Session) 


BEGAN RECORD ENG SORSIERMINAL SESSION 


R; T=0. 


O170 502.6330 223 


graphics 
EXECUTION BEGINS... 


DO YOU 
1. 
Ze 
36 
4. 


OPTGRAPH 


DESIRE DOANE. b= 
OPEN LOCP TRANSFER FUNCTION 
NOISE TRANSFER FUNCTION 
COMPENSATOR TRANSFER FUNCTION 
PATA OF TGRAPE 


ENTER OPTION NUMBER. 


3 


OPTGRAPH 


46 


THE COMPENSATOR TRANSFER SSY STROOD] 2 Dae 
OPTSYS CChwaeais: 


le. NUMBER SOF “STATES ]— Z 
Za NUMBER OF CONTROLS (1NPUtS— 1 
3. NUMBER OF MEASUREMENTS S(CUIC Ut ae 


4. MARKOV PARAMETER = 10.0**- 6 


CLEAR SCREEN TO COUT 


OPTGRAPH 


AN OPEN LCOP POLE-ZERO, ROOT LOCUS, BODE, NYQUZST, 
AND/OR NICHOLS PLOT IS DESPier- fon. 
INPUT # ? 


CUTPUT # ? 


OPTGRAPH 


COMPENSAICR TRANSF CR FUNCT Ben 
INPUT # = 1 
OUTPUT t= = 
DO YOU DESIRE GRAPHICAL RESPONSE AND/OR TABULAR 
DATAPEORG 
1. FOLE~ZERO MAP 
Ze ~ROOT-ECGCIS 
3.) SEGUE 
4. NYQUIST (POLAR PLOT) 
o- NICHORS 
6. ANOTHER INPUT/OUTPUT COMBINATION (SAME 
TRANSFER FUNCTION) 
7. ANOTHER DYPE TRANSFER RUNG? LON e7 eee 
CPTGRAPH 
ENTER OPTION NUMSEEe 


4°] 


BODE FLO. 

SCOOP ENSA Terror fhe tJ NCL ON 
DOT VOU Prosi kks 

1. GRAPHICAL OUTPUT 

2. TABULAR DATA 

Sa CULTS UBROULINE 

ENTER OPTICN NUMBER. 
1 


BODE PLOT 
COMPENSATICR TRANSFER FUNCTION 
PLOTTER SELECTION AND PAGE SIZE 
(NOTE: PAGE IS 11.0 INCHES BY 8.5 INCHES WITH 
SELECTICN OF TEK618) 
WHICH PLOTTER DO YOU DESIRE: 
1. TEK618 
2. VERSATEC 
CHOCSE OPTION 1 OR 2 


2 
PAGE SIZE (MAXIMUM = 21.0 INCHES BY 21.0 
INCHES) 
HEIGHT = 
5.65 
WIDTH = 
S255 


PRO EIeNG LIN ETSOPORSGRAPH 
FREQUENCY RANGE - RADIANS PER SECOND 
MINIMUM FREQUENCY 2 


0.1 


Ui 
“) 


MAXIMUM FREQUENCY 
1000000 


PLOTELENG LIMits FOR GRAPH 
PRE@UENCY RANNGE = "RADIANS PER SECONDE 
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MEN Utes OU Ne 

MAXIMUM FREQUENCY 

DO YOU DESTRE TOMA KE ANY eG HANG ie 
TYPE "Yes OR aia. 


a 


PLOT HEADING 


0. 100E+00 
0. 100E+907 


HOW MANYoL ENS OF SBEADING DOV Oy Distr 


( 3 MAX ) 
ENTER NUMBER OF LINES. 


{ ENTER “OUSFOR NOr Hea EiliG y) 


Z 
A MAXIMUM OF 32 CHARRC LER S PER baie 
ALLCWED 
LINE Ios 
filter sinrulation 
DN Zee: 


ref. 8; ~ages 352° =s524 


PLOT HEADING 
LINE tes FILTER SIMULATION 
i ENE REF mes PAGES 332 - 
DO YOU DESIRE “ORMAKE ANY CHANGES) 2 
TYPE "YES" OR “Nege 


334 


A DISSPLA METAPILE HAS Bary Crates 


CLEAR SCREEN TO COM: 


GAIN MARGIN AND PHASE MARGIN 
PHASE CROSSOVERS ERE UeNGy oS 75.1184 
GAIN MARGIN I5 =5 9 SOO DE 
GAIN CROSSOVERS PREOQULNCY [Ste 
PHASE MARGIN IS -89.9480 DEG 
CLEAR SCREEN TO CONTINUE 
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RAD/SEC 


KAD/ 


" 
tad 


C 
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BODE PLOT 
COMPENSATICR TRANSFER FUNCTION 
DO YOU DESIRE TO MAKE ANY CHANGES TO: 
1. ELOTTER / PAGE SIZE 
2. GRAPH LIMITS (FREQUENCY RANGE) 
3. HEADING 
4. NO ADDITIONAL CHANGES - PLOT BODE PLOTS 
5. NO CHANGES — EXIT BODE SUBROUTINE 
ENTER OPTICN NUMBER. 
s) 


BODE PECL 
COMPENSATCR TRANSFER FUNCTION 
PCy ,OUSDESE RE: 

LeohnnenicAL CUTPUT 
2- TABULAR DATA 
Seoul. SUBTOULLN & 
ENTER OPTICN NUMBER. 
3 


OPTGRAPH 


COMPENSATCR TRANSFER FUNCTION 

INPUT # 1 

OUTPUT # 1 
DO YOU DESIRE GRAPHICAL RESPONSE AND/OR TABULAR 
DATA FOR: 

1. EOLE-ZERO MAP 

2. ROOT-LOCUS 

3. EODE 

4. NYQUIST (POLAR PLOT) 

5. NICHOLS 

6. ANOTHER INPUT/OUTPUT COMBINATION (SAME 
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TRANSFER FUNCTION) 
7. ANOTHER TYPE TRANSFER FUNCTION /EXIT 


CPTGRAPH 
ENTER OPTION NUMBER. 
Fi 
OPTGRAPH 


DO YOU SDESITRE TOmalaiyZ E : 
1. OPEN LOOP TRANSFER FUNCTION 
Ze NOISE TRANSFER FUNCTION 
3. COMPENSATOR TRANSFERS FUNGria 
q. EXIT OPTGRAPH 

ENTER OPTION NUMBER. 

4 

ReOT=720 176 54.015: 39 s00 

LeCOEd OLE 

END RECORDTING OF TERMINAL SESSION 


The preceding example follows at the end of the 


chapter as figure 3.3 and figure 3.4 . 
b. Example of Tabular Data Listing File 


The terminal session for the Bode tabular data response is 
Similar to the Bode graphic reSponse terminal seSsion. 
Appendix D contains an example of tabular output sent toa 
disk (OPGRAPH LISTING) for five points between 10 and 100 


radians/second. 


3. Nyguist Respcense 


The graphic response section of the subroutine NYCST 
creates a polar plot from the magnitude and phase data 
received from the sukroutine FREQ. The plot 2s auconma ti-— 


cally scaled from a minimum radius of one by using the 


Si 


Subroutine MINMAX to determine the maximum phase and magni- 
tude values. The gain / phase margin calculations and the 


tabular data output are identical to the subroutine BODE. 


4. Nyquist Respcnse (Interactive Example) 


The following example is of the open loop transfer 
function for the filter simulation example described in the 
Bode interactive section. 

This example follows at the end of the chapter as 


mrgure 3.5 . 


5. Nichols Respcnse 


The graphic response section of the subroutine 
NICHCL creates a single rectangular plot of magnitude in 
decibels versus phase in degrees. The magnitude from the 
subroutine FREQ is converted to decibels and plotted with 
the corresponding magnitude. The plot is automatically 
scaled using the subroutine MINMAX to determine the maximun 
phase and magnitude values. The gain / phase margin calcu- 
lations and the tabular data output are identical to the 


SubLoutine BODE. 


6. Nichols Response (Interactive Example) 


= == 2 SSeS ae ee cee ee ee = =a oS == 


The following example is of the open loop transfer 
function for the filter simulation example described in the 
Bode interactive section. 

This example follows at the end of the chapter as 


Prgure 3.6 . 
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IV. CONCLUSIONS AND RECOMMENDATIONS 
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A. CCNCLOUSIONS 


An evaluation of the computational ability of OPTISYSX 
and CPTGRAPH, the frogram was tested using an 82 X 82 
matrix, provided by NASA-Edwards, of longitudinal eguations 
for the X-29A experimental forward-swept wing fighter 
ale CcCrait pPLotoby pe. The OPTSYSX program array had to be 
redimensioned anda 2-megabyte virtual machine size was 
required for a system of this size. 


The open loop eigenvalues calculated by the OPTSYSX 


program compared favorably with the Eigensystem data 
supplied by NASA-Edwards, but the OPTSYSX program failed in 
calculating the open roo transfer function ZELOS. 


Extraneous zeros (calculated by OPTSYSX) were of the same 
approximate magnitude as the open loop system Zeros. This 
prevented the identification and elimination of the extra- 
neous zeros by comparison with the Markov parameter. The 
open loop zero identification was further complicated by not 
having available open loop transfer function zero informa- 
tion from NASA-Edwards to compare the OPTSYSX results to. 

When reguested for additional information NASA-Edwards 
frrovided a computer tape with a revised 98 X 98 matrix of 
lon g1-t0 dana Laio< Ion equations for the X-29A fighter 
alceratt. For this revised system, the OPTSYSX progran 
requires a virtual machine size in excess of 2-megatyte. 
Time constraints have prevented running the OPTSYSX program 
with the revised systen. 

The OPTGRAPH program was tested using the OPTSYSX 
Frogram's open loop Eigenvalue and transfer function cutput 
data for the X-29A 82 XK 82. Systen ct Wongitddinat nortan 


One, 


equations. OPTGRAPH provided excellent quality Bode flcts 
for all of the open loop transfer function inputs and 
outputs and did not exhibit any computational difficulties 
as indicated by the lack of error messages at the end of the 
test. While the test results are encouraging as to the 
ability of the OPTGRAFH program to assist in the analysis of 
large order systems, the actual computational accuracy of 
the OFTGRAPH program can not be fully verified until accu- 
rate transfer functicn information can be obtained from the 
Gra srSk progran. 


B. RECCHMENDATIONS 


Based on the results of this Thesis the following areas 


emerged as possible areas for further research and study. 


1. Root Finding Routines 


The OPTSYSX frogram has two similar root finding 
Subroutines (HQR and HQR2), to find the transfer function 
EFlgensystem (poles) and zeros. The order of the transfer 
function numerator is calculated separately in the sutfrou- 
tine ZEROS. 

The inclusion of the open loop transfer function 
elgensystem (pole locations), zero locations, transfer func- 
tion numerator orders, and transfer function gains, for the 
revised (98 X 98) lcngitudinal mode system by NASA-Edwards 
should simplify the lcecating the problem in calculating the 
transfer function zeros. 

The open loop system data from NASA-Edwards also 
included Bode magnitude and phase plots, and gain and fhase 
Margin information, which can be used to evaluate more 
conclusively the OPTGRAPH program's ability to help evaluate 


a large crder system. 
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2. Program Memory Requirements 


> a = == [= a eee Ss 


The OPTSYSX frogram reguires in excess of 2-megabyte 
of virtual machine memory when configured for large matrix 
operations (98 xX 98). Virtual machines with this remory 
capacity are not normally available toa user. The memory 
usage for matrix storage and manipulaticn 1s a possible area 
for the reduction in program memory requirement size. The 
OPTSYSX program retains most of its computational arrays in 
memory. This method simplifies programming but is extremely 


inefficient in its use of virtual machine memory. 


3. Further Modifications 


= — a EE ME ee ES Se Se ee Se 


The area of modern digital controls should be further inves- 
tigated. The computational abilities of the OPTGRAPH 
program make it readily adaptable, with minor modifications, 
for assisting in the analysis of transfer functions in the 


digital domains. 
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GEPSYs GAEC LISTING 


Een ACE OFF 
ea kk kok tok kkk kok kk kk do Rk kok kk ok ok kok kk ok 
Dimmer Lot AL © 
SONOS POLL SYSX,, ORPTCALC ANDSOPTPLOT 
Peo wince lie ME RESPONSE OFFA SYSTEM. 
Depts ae DL EL 
Vhs EON 1.0 16 JULY 1984 
POUT TED Sri testes uae TAS 
He HK RR HB RK oo ok ok oko oo Rk oo ok ok 
SiROoOCEN 
worG@ixet -~ENDFOUR 


+ «© + + + & 


He He He Ee he He ie oe ee oe eK ee 2 ee eo 2c 2h fe oe oe fe oie ee 2c ee ie oe a 2c i a i ee ie oe ek ok 


CUtpde FROMMOPTSYSX WiLL COME TO THE SCREEN IF YOU 


WLSH ELIS K. PELE 
ENTER OD PNY SOL ke eet hey LELD SCREEN 


He He 2 ke ke fe ok oe ee he ee ee 2c ee a ee 2 ee oe ee eo he ee he eee ae 2c ic ee a ee eK oe Ke ee a ee 


-ENDFOUR 

GSREAD VAES &ANS 

SIF .&ANS NE .D &GOTC -ONE 

PIGEDEF 06 DESK OUTPUIX LISTING Atl 
CLRSCEN 

EBEGTYPE -ENDSI2X 
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FRR RRR RRR RR Ra RR RRR RO RO a ROR RO RR Rk ka ka Rk ke at aka ica 


OUTPUT WILL GO TO DISK FILE *OUTPU DYES eee 


We Hee ee He he he ee he ee ee ee ee ae ee ee ee eo 


=ENDS DX 

CP SEEER +s SEG 

=ONE 

CLRSCRN 

GCP L&aNK O0039P Wale Zo ieee 

AGGeZEtT F 

FOR RO CR RR REI IIR RR IRR IRR IR RR RR Ra aR aC Rake aie a ote a ake Fe ak 
* * 
* ALLCW THE USE CF AN OLD “OPTMAT DATA A1" “ 
* * 


SOCIO III II GI IOITOI IGG IG RII Ia a ii a kak kok kok a aOR ok ake aie okeak ak ae 
RENAME OPTMAT DATA Al OPTSYS DATA Al 

SIF GRE NE O &S&GOTO =FIRST 

RENAME OFTSYS DATA Al OPTMAT DATA Al 

FINDSTAK OPTMAT DATA A1C0O1 0 LIMO002 ALL GROUP1 

SREAD VAKS &A1 &A2 

SIF .€A2 NE .1 &GOTO -START 

CLRSCEN 

SBEGIY PE. — fine a we 


We MMe Re He Re Ee ee ee a a Re ee ee ee ee oe a ee oe eae a eo eo a 


YOU HAVE A DATA FILE NAMED “*OPTMAT DADA ey oune: 
DISK THAT WAS PREVIOUSLY GENERATED Bieges Orton. 
FROGRAM AND CCNTAINS THE F, G, H, GAMMA, A AND B 
MATRICES FROM THAT RUN. 
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Taerou fwOUDD LIKE STO SUSE THESE SAME MATRICES FOR 
oer [he CPISYS PROGRAM “WILL ) READIN. THE 
DESIRED DATA AT THE APPROPRIATE TIME, 


IF YOU TYPE (Y) ES. 


Peo cnnke NP UT Moe EouUET IN THAT FILE BEING ERASED! 


eA A Re aR i ek aR aR a a a Re ok kee ke a ae a eo ae ook ok ok aK oe ko ok ok ok kok oko ok ake ok 
-ENDTWO | 
SREAD VARS &ANS 

PeeeoANs EO .¥ SGOTO —START 

eR OL 


took kok ok ok doko kkk kkk ok doi k kok ok kkk kkk kok kkk dk kkk kok Ok ok 


FRASE THE OLD “CPTMAT DATA Ai" DATA FILE 
peer O00 O0' VIN THE NEW “ORTMAT DATA FILE" 
TO ACT AS A FLAG FOR OPTSYSX AND OPTCALC 


+ + + &+ + 
+ + + + + 


ee ee Re RR KR ke Ce a a RR a RR oR oo oe a a RR a a Ro ke a aK oe kek ook 
ion OPIMAT DATA AT 

GCaencK 000 0 

iiemortcK OPIMAT DATA Al F 80 1 

SOLA RT 

CLRSCEN 

SBEGTYPE -ENDONE 


mabeOrlsyYS EXEC GONTROLS B@URS PROGRAMS: 


1 OPTSYSX FCRTRAN (SYSTEM ANALYSIS) 
GENERATES OPTMAT AND OPTGRAPH DATA SETS 
= RERUN USES OPTNAT 

POE TORE GMLCHIRAN (CALCULATE IINE RESPONSE) 
GENERATES OPTPLOT DATA SET USES OPIMAT FROM 
OPiS 1 S:X4 

3 OPTGRAPH FORTRAN {(POLE=-ZERO, ROOT-LOCUS, BODE, 
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NYQUIST, ETC) 
USES OPIGROL -OPTGRNO -OPTGRCM FROM OPTSYSX 
4 EXIT 


OPTPLOT FORTRAN 2S _ THE FOURTH Se ReGnA 

USES OFIPLOT DATA SET TRON CrPreme 
THE SIZE OF THE DATA Stet S VA RY oe eee 
ORDER, AND CAN USE ABOUT 20% OF THESUSERSSDLSKSPAGre 
THEREFORE ENSURE THAT SUFFICIENT DISK see aeteer. 


AVAILABLE. 
IF DATA SET IS ALREADY AVAILABLE YOU MAY RUN ANY 
PROGRAM 
- ENTER 1,2, 3, 4 ANY OTHER INPUT 
RETURNS TO MENU 
-EN DONE 


SREAD VAKS &ANS 

SIF .6ANS EQ 21 6GOTO] Ob Tas 

CIF .GANS EC 22 €CGOT0 —erElerre, 

SIF .&ANS EQ .3 &GOTO -OPGRAPHI1 

SIF .~.&ANS EQ .4 GEEXIT ERC 

&GOTO -START 

=OP is : 

FILEDEF 8 DISK OPIPLC] DATA Alea 

FILEDEF 9 DISKSOPIMAT DATA 21a 

FILEDEF 10 DISK OPTGRCL DATA Al (PERM 

FILEDEF 01 DISK OPTGRNO DATA Al (PERM 

FILEDEF O04 DISK OPTGRCM DATA Al (PERM 

GLOBAL TXATILIB VFORTLIE INSEDe NCueysE 

STYPE LCADING OPTSYS.... GENERAL) SP Ad bveseet oclet ens eae pone 
PROGRAM 

STYPE.... OPTMAT DATA Al 2.246. INPUT Sta sol ON Oe 

OPT S We 

CLRSCRY 

STYPE... DATA SEIS “GeeGrtaAt AND .-. OPGRAPH CREATED 

CP SPREE Rs aoe 
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wor ~SIART 
~OPTCALC1 


FE FI He RT Re ee ee oe oe ie ae ie ie ae I a ote ae oe a ae ae ae fe ke ake ae fe ok eae oc oe ee ok ae ate ote ote ke ote ok ok ok ok ake ake ake 


x x 
* CHECK FOR DATA IN THE FILE "OPTMAT DATA " a 
* BEFORE LOADING OPTCALC * 
* x 


AE HE Ke AE Ie A eR a 2 ae ke ik ae aie ee ae tee a oe afk ae fe fe a ae ik a aie fe fe kc fe eo kc fe ae oie fe fe oe ake oe ak kc fe ee fe ae 
FINDSTAK OPTMAT DATA Ai CO1 0 LIMOO2 ALL GROUPT 
SREAD VARS &A1 &A2 

SIF .€A2 EQ .1 &GOTO ~-OPTCALC 

STYPE PROPER DATA FILE IS NOT AVAILABLE FOR OPTCALC 
SPSSLEEPS> SEC 

ecemo —SIART 

=OPm@ CALC 

PPM@E oiaisis oo so OP TCALC.... TIME RESPONSE PROGRAM 
PiGeosr "2 DISK OPTPLCT DATA Al (PERM 

poorer 9 DISK OPIMNAT DATA Ail (PERM 

SUObrmrE LXTLIB VFORTLIE IMSLDP NONIMSL 


CELTEALC 

STYPE -- OPTPLOT CATA Al CREATED 
Bemoreee 5 SEC 

CiRSCRHN 


SBEGTYPE -ENDNINE 


He HK ee ae ae ate fee ae ae oe aR 2 oe 2c ee ae ok ae fe ee oie ie fe ae aie i i oie ae eee ae ie ok ke ok ae ice ac ae Ae oe ae ie ae Ie 3K 


me YOU ARE SMSSHTISFILED WITH THE RESULTS 
PAUSSEAR AND WOULD” LLKE TO EXIT TO Cus, 


eee eer LO EX Lis 
(ANY OTHER INPUT TO CONTINUE) 


66 


He Ko eK a he A A he KK i i eK fe he a KK A KK KK KK HK KK 


~ENDNINE 

GSREAD VARS 6&ANS 

SIF .&ANS EC . YRCCOT GR 

crYPE ELOTSANe wa 

CP SLEEP 4 SEC 

GLOBAL TXTLIB DISLIBVS S2ZDISLIB S2ZINTLEB VT ORE ee Gre Lie 
NONDTESE Cis is 

FILEDEF 8 DISK OPTPLOT DATA Al) (OSE 

FILEDEF 11 DISK LBLANK IMP AS (RECFM F LRE@E 2000 5p a 2 
2400 PBLKSIZE 2400 XTENT 600 

FILEDEF 12 DISK EISSCCMAP MAPDTA * (REGEM Perkeereccs 
BLKSIZE 400 

FILEDEF 13 DISK TABLET IMP A3 (RECFM VS LRECL 2038 BiWsiZ5 
208 XTENT 1000 

FILEDEF 14 DISK LBLANK MAPDTA SLBMODE (RECFM VS LRECL 608 
aE SARA, UIs: 

FILEDEF 17 DISK DISSFIAS@SYSUL! AS eGR ECE Sewer Bere) se 
BLCGK. 2000 XTENT aod 

FILEDEF 18 DISK DISSEIA METAFILE AG (RECPM VESSLRECL 120c8 
BLEEK I069 


OPTPLOT 

&GOTO -META 

-OPGRAPH1 

eK He He He He He He he he He a he eK He ee he KK oe eK kK eK KK ok KK KK KK 
‘fe ae 
2 CHECK FOR DATA IN THE OPTGRAPH DATA SETS " x 
* EEFORE LOADING OPGRAPH = 
® x 


Ae SK aK eee a a a RR a a a a a oo oR oa ao oa oe kok oe ck ek 
2 eke ee a ee a oe a oa eo oo oo ok tk eK 
RENAME OFTGROL DATA Al OPTSYS DATA Al 

SIF GRC NE 0 €GOTO —CEiNoms: 

RENAME OPTSYS DATA Al OPTGROL DATA At 
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FINDSTAK OPTGROL DATA A1CO1 0 LIM002 ALL GROUP1 

SREAD VARS &A1 &A2 SA3 SAY SAS 

SIF .&A2 NE .0 &GOTO -SIXTEEN 

-OPINCISE 

RENAME OFTGRNO DATA A1 OPTSYS DATA A1 

SIF &RC NE 0 &GOTO -CETCOMP 

RENAME OPTSYS DATA A1 OPTGRNO DATA A1 

FINDSTAK OPTGRNO DATA A1 CO1 0 LIMO02 ALL GROUP1 

SREAD VAES 6A1 &SA2 &A3 SA4 EAS 

SIF .6A2 NE .0 &GOTO -SIXTEEN 

-OPTCCMP 

RENAME OPTGRCM DATA Al OPTSYS DATA Al 

SIF &RC KE 0 &GOTO -NCDATA 

RENAME OPTSYS DATA A1 OPTGRCM DATA Al 

FINDSTAK OPTGRCM DATA A1 CO1 0 LIMOO2 ALL GROUP1 

SREAD VARS 6A1 &A2 &A3 SA4 EAS 

SIF .&§A2 NE .0 &GOTO -SIXTEEN 

-NODATA 

SPRINT PROPER DATA FILE IS NOT AVAILABLE FOR OPGRAPH 

CP SLEEP 5 SEC 

EGOTO -START 

-SIXTEFEN 

STYPE ...OPTGRAPH.... CLASSICAL ANALYSIS OF OPTSYS OUTPUT 

SET LDRTELS 10 

GLOBAL TITLIB DISLIBVS 92DISLIB 92INTLIB VFORTLIB GRFLIB 
NONIMSL CMSLIB 

SET LDRTBLS 10 

FILEDEF 01 DISK OPTGRNO DATA A1 (PERM 

FILEDEF 02 DISK OPGRAFH LISTING A1 

FILEDEF 03 PRINTER (RECFM FA BLKSIZE 133 PERM 

FILEDEF 04 DISK OPTGFCM DATA A1 (PERM 

FILEDEF 06 TERM (RECFM FA BLKSIZE 133 

FILEDEF 05 TERM (RECFM FA BLKSIZE 80 

FILEDEF 10 DISK OPTGRCL DATA Al (PERM 

FILEDEF 11 DISK LBLANK TMP A3 (RECFM F LRECL 2400 BLKSIZE 
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2400 XTENT 600 

FILEDEF 12 DISK ISSCCHAP MAPDTA * (REGEN fF Ur oehwea 
BLKSIZE 400 

FILEDEF 13 DISK TABLET TMP A3 (RECFM VS LRECL 2ZOGe ri: 
20G) XIENT 1008 

FILEDEF 14 DISK LBLANK MAPDTA ELBMODE (RECHUBVSG SLE Eee. 
Bik ool Ae S 

FILEDEF 17 DISK DESSFLA SYSUT1 HS (RECPRETE Samael 
BLOCK 2000) YSN T 500 

FILEDEF 18 DISK DISSFIA METAFIEE AY (RECrM Ves CRP ere Uc 
BLCCK 19069 


CPTGRAPH 

-META 

ete ee RCRD ea fe ok EA ok shen doko oe deevck ok ac ete cone see petal aes kore 
% = 
® CHECK FOR FILE "DISSPLA METAFILE A4" ON * 
. THE USER'S DISK BEFORE GOING TO DISSPOP * 
* =~ 


FOC GIORGI ok a a a kok ak a ok ak ok ak ake okie ok ak ak kc ofc ok kk ofc ak otek ic ak ofc ok dk ak ke ok ak ake ke ak 
RENAME DISSPLA METAFILE AY OPTSYS METAFILE A4 

SIF &RETCODE NE O &GCTO -START 

RENAME OPISYS METAFILE AY DISSPLA METAFILE A& 

=T Gin 

CERSCEN 

EBEGTYPE =ENDTEN 


Ae AK AK ee i ei ee ee I ek ok eK i oe ek a ke ee kK OK ee IK Ke eK KK KK KK 


DC YOU WANT A VRSTEC PLOTTER SHOOTH Cen eCre ta: 
THE DISSPLA METAPSLE THAT YOU JUS tere eee 
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(Y OR N) 


ok kK a kok ak gk ok kkk ko kkk kok kook kk koko ak ak tok took kok ook 


=ENDIEN 

SREAD VAES &ANS 

SIF .&ANS EQ .Y &GOTC -NINTH 
GifmeechNs EO .N &GGOTC -START 
Geom —EIGHT 

=Ni NTH 

Eee OLSSPOP VRSTEC 


Cin ChN 
CeeGiYPE —ENDTWELVE 


THK IR IR ee kK 2 KK A OK a OK Ae ak a ok RR ok oe ak oe ke ik ak ie 2k 3K a ok ak aK ak ae A oie OK a a OK ik ie ok oe ok a eK kak ok ok 


peu GRaArH(S) CAN BE PICKED UP AT THE COMPUTER CENIER. 
Miterciumeen (>) WElLl BE ADDRESSED TO “POP (USER ID)". 


EK KK OK OR eo ee aK OK oe a ie oie ie ie eK ie Ok ek OK ie oe ie a ie ek ak ie ak Ie OK ie SK OK OK OK eK OK eK ok Kk oe ok ke kK KK 


~ENDTWELVE 
Smo LE bP 5 SEC 
SGOl0 -START 
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APPENDIX B 
OPTSYSX PROGRAM LISTING 


This portion of the thesis contains the OPTSYSX FORTRAN 
program (88 pages). 
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APPENDIX D 
OPGRAPH LISTING 


This portion of the thesis contains a sample of tatular 
output sent to a disk (OPGRAPH LISTING) 
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